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Abstract 

We review the theoretical model and experimental realization of the atom optics 5— kicked rotor 
(AOKR), a paradigm of classical and quantum chaos. We have performed a number of experiments 
with an all-optical Bose-Einstein condensate (BEC) in a periodic standing wave potential in an 
AOKR system. We discuss results of the investigation of the phenomena of quantum resonances 
in the AOKR. An interesting feature of the momentum distribution of the atoms obtained as a 
\ result of short pulses of light, is the variance of the momentum distribution or the kinetic energy 

ijp') /2m in units of the recoil energy E^ec = hujrec- The energy of the system is examined as a 
function of pulse period for a range of kicks that allow the observation of quantum resonances. 
In particular we study the behavior of these resonances for a large number of kicks. Higher order 
quantum resonant effects corresponding to the fractional Talbot time of (l/4)Tj' and {1/5)T't for 
five and ten kicks have been observed. Moreover, we describe the effect of the initial momentum 
of the atoms on quantum resonances in the AOKR. 
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I. INTRODUCTION 



In classical dynamics, small perturbations in the initial conditions of a system may grow 
exponentially with time and produce large ones in the final states. This makes the motion 
practically unpredictable, and the phenomenon is known as classical chaos On the 

other hand, the investigation of the behavior of quantum systems, the classical limit of which 
are chaotic, has brought up a new discipline in physics, which is now known as "quantum 
chaos". The pioneering work by Giulio Casati and co-workers predicted that a quantum 
particle shows diffusion following the classical evolution, after a certain time known as the 
quantum break time . The diffusion stops due to quantum interference, beyond the 

quantum break time, leading to dynamical localization js], Q] which is analogous to Anderson 
localization in solid state physics. 

The strange connection between classical chaos and quantum chaos has forced people to 
develop models that can describe systems that are chaotic. The kicked rotor is one of the 
simplest models that is useful to understand classical and quantum chaos. The classical 
model of the kicked rotor is a simple pendulum that is kicked periodically by an external 
force. The dynamics of the the pendulum that exhibits chaos can be described by the 
Hamiltonian mechanics. Classically a kicked rotor is a well known system that makes its 
quantum analogue the best model to study quantum chaos [s ill. 

The quantum analogue of the classical kicked rotor is obtained by a cloud of ultra cold 
atoms interacting with a pulsed standing wave of laser light and is known as atom-optics 
kicked rotor (AOKR). The atom-optics realization of the kicked rotor was first demonstrated 
by the Raizen group js] in 1995, and has been studied extensively theoretically and experi- 



mentally in 



18| and 19N28| and the references therein. The AOKR model has allowed 

quantum resonances" 



experimental studies of the phenomena o: 



the matter-wave Talbot effect 
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29| that occur due to 



30H32|], analogous to the optical Talbot effect |33|. 
An interesting feature of the momentum distribution of the atoms obtained as a result 
of short pulses of light, is the variance of the momentum distribution or the kinetic energy 
(p^)/2m, in units of recoil energy Erec = ^'^rec = 'h?'k\/2m, where = 27r/A, where A i 



IS 



the wavelength of the laser beam and m the mass of the atom. The energy of the system is 
examined as a function of pulse period. 

In this paper, we present some of the known results along with our new results to illus- 
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trate the phenomena of quantum resonances in the AOKR. We also present experimental 
observation of higher order quantum resonances corresponding to the fractional Talbot time 
: br five and ten kicks. The results are in good agreement with those presented in Ryu et. al. 



30| . where higher order resonances were found in a similar experiment with a condensate of 
Na atoms. Moreover, in section III, the dependence of the initial momentum of the atoms 
on quantum resonances have been discussed in detail. The results presented are in good 
agreement with those proposed in [l^ . 

II. THEORY 

The general Hamiltonian that describes the dynamics of atoms in the Atom-Optics kicked 
rotor experiment can be expressed as 

2 ^ 

H = ^ + n<i)dCos(2kLx)^5{t-nT), (1) 

n=l 

where T is the kick period and 0^ is the kick strength. The kick strength 0^ = ficffT, which 
is determined by the effective Rabi frequency during the laser pulse VL^^ and the pulse length 
r. The effective Rabi frequency is defined as f2efr = ^^^/A. 

A. Evolution of the wave packet 

The time evolution operator that determines the evolution of the wave function from one 
kick to immediately after the next kick is known as the Floquet operator J-" {27! . Issl . In 
general, the state of the system \ip{tn)) after n number of kicks can be found by the Floquet 
operator acting on the initial state \ip{tQ)) i.e., 

m^)) = :F-m,)) . (2) 

In the first part of the Hamiltonian describing the quantum kicked rotor, the system evolves 
freely via the operator 

Ufree = e'^^'^/'^\ (3) 

Then delta kicks act for a small period of time, and the evolution is given by the operator 



3 



Combining the above two operators, the Floquet operator is given by 

■p _ ^-ip'^T/2mh^-i(j)acos{2kLx)T/h ^g-^ 

^ free ^kick 

or in a more general form as 

J-" = UfreeUkick- (6) 

From the definition above, the action of each kick on the atoms, which is governed by Ukick 
is followed by a phase of free evolution Ufree- A more detailed description of the evolution 
of the quantum wave packet for a delta kicked rotor system can be found in jl^ . 

B. Simulating the AOKR 

In order to simulate the AOKR, the well-known split operator method for the time 
evolution of the Floquet operator has been used, which is ideally suitable for the experiment 



with 2^^ grid points, as described in detail in 27|] and 28[. As described in the previous 
section, the time evolution of the quantum 5— kicked rotor can be separated into two parts: 
The evolution of the system during the time when laser pulses are applied on the atoms as 
kicks, and the evolution during the time of free expansion in between the kicks. The kick 
potential is a diagonal operator in position space, whereas the free evolution is a diagonal 
operator in momentum space. Therefore, for the time evolution of the system, a Fourier 
transformation is used in order to transform the wave function from position space into 
momentum space and back again. The initial wave function representing the state of the 
system can be defined as 

The width of the initial wave packet is chosen such that it covers several wavelengths of the 
kick laser. 

For the kicking part of the time evolution, the delta-function is represented by a block 
function with a width r <^ T and an area 0^ |27|. The initial wave function tp{x) in position 
space in the simulation is multiplied by 

exp{^}. (8) 

r d 

The wave function is then transformed from position space into momentum space repre- 
sented by ^{x) using an inverse Fourier transform. The free evolution part of the system 





Pn (recoils) 



FIG. 1: The simulated momentum distribution of the wave packet for the quantum resonance case, 
i.e, ^ = 47r is illustrated for each kick, with an initial p = momentum state. The distributions 
are shown for = 1 — 6 kicks, with a kick strength of (pd ~ 1-5. 



can be determined by multiplying the wave function in momentum space by 



(9) 



which gives the final state in the momentum space (p{x). The parameter Tf^ee represents 
the pulse period in between the kicks. After the free evolution part, the wave function is 
transformed back to position space by applying the Fourier transform and gives the updated 
state tplx). The process is repeated for each kick, giving us the updated wave function each 
time after the evolution. The momentum and energy of the atoms after each kick is then 
determined. 
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C. Quantum resonances in AOKR 



Quantum resonances are well known quantum effects that occur for specific initial con- 
ditions and certain values of the kicking period. We investigate these resonances in the 
AOKR by considering the atoms diffracted by a standing wave potential. A "quantum res- 
onance" exists for T = rriTT, Tt being the Talbot time, where all kicks add coherently and 
energy grows quadratically with the number of kicks. An "anti-resonance" is observed at 
T = (m + |)'7r/2wrec5 where the effect of each kick is effectively negated by the following 
kick. 

The free evolution term is of primary importance in describing the quantum kicked rotor. 
If the initial state is assumed to be a momentum eigenstate |m), and let the operator Ufree 
act on it such that /)|m) = mklrn), then the evolution is different for different values of ^. 
If k is an even multiple of 2tt, then the free evolution results in unity. This means that the 
free time evolution in this case does not have any effect on the state vector of the system. 
The system responds only to the 5— kicks js^, and the evolution is fully governed du ring 



interaction with the kicking potential. This is the condition for quantum resonance |12 |. 
On the other hand, if k is an odd multiple of 27r, then the initial momentum eigenstate 
acquires different quantum phases for even and odd m. The state after the free evolution 
acquires a quantum phase of +1 for even values of m, whereas for odd m the quantum 
phase results in —1. For odd values of the momentum components of the wave-function, the 



herefore, ^ = 27r is known 
SgI . [3^. The phenomenon 



phase of —1 leads to an oscillation of the energy of the system, 
as the quantum anti-resonance condition for the kicked rotor 
of quantum resonance and quantum anti-resonance is demonstrated in the simulations in 
Figures [1] and [2] respectively. 

The general expression for the qth order momentum moment of the 5— kicked particle 
after time nT as expressed in [U] is given by 



j=-oc 

The parameter T is defined as 



sin(?2T) ^ "' ^ 
sm I 



(j + pr- (10) 



T = ivr(l + 2/3)/, (11) 

where /3 is the quasimomentum. Quantum resonant and antiresonant behavior of the system 
can be characterized by the second order momentum moment for g = 2, and is proportional 
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FIG. 2: The simulated momentum distribution of the wave packet after kicks for the quantum 
anti-resonance case, i.e., = 27r is shown for each kick. The parameters are = 1 — 6, with a kick 
strength of w 1.5. 

to the kinetic energy. 

As discussed in the previous sections, the evolution of the system can be determined by 
the time evolution operator. In the case of quantum resonance, the free evolution part is 
unity, therefore for kicks the evolution operator is written as 



u 



(12) 



where n/Ti = (p^ is the kick strength. The final state of the system is as if a single 5— kick 
had been applied with a strength times the original kick strength 0^. The probability for 
an initial momentum state |m), to be populated in the final state \n) after kicks is given 
by 

(13) 
(14) 



n\U^\m 
n\e^ s ^ >)\m 
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Using the definition of orthogonality of the momentum eigenstate {n\m) = Snm and the 
identity e^*'^^™^') = J^e*''^, where J is the ordinary Bessel function of order s, we get 

Pn = Jn-m (^~^^ ' ^^^^ 

For a system that starts in a zero momentum eigenstate, the final momentum distribution 
becomes 




The energy of the ensemble after kicks can also be found at the quantum resonance, and 
is given by 

En-Y.^^'Pn-~. (17) 

n 

For a kicked rotor starting in the zero momentum state, the above result shows a quadratic 
growth in energy for kicks for the quantum resonance. 



D. Experimental effects 

In the experiment, the initial momentum distribution from the BEC has a width of a ~ 
0.18 recoils. The initial width associated with the BEC qualitatively alters the momentum 
distribution after the kick sequence. In an ideal case the initial momentum distribution is 
a delta function at pj, however in reality, the initial momentum distribution has a width. 
Because of this width, there are some atoms at a momentum pj + Ap, where the quantity Ap 
is much smaller than one recoil. In order to take into account this initial momentum width, 
we performed the simulation for a range of initial momenta. In Fig. [Ill the probability of 
finding the atom at a final momentum is plotted for a range of initial momenta for both 
antiresonance (top) and resonance (bottom). The initial momentum distribution of the BEC 
is also displayed (as a solid curve). The effect of the initial momentum on the system can 
be observed as a transfer of the probability of atoms to states at small offset momenta. As 
the probability is very sensitive to the initial momentum, this sensitivity therefore translates 
into peaks appearing at a small offset momentum pi + Ap which are not present at pj, where 
Pi is equal to one recoil in this case. Conversely, this can also be seen as a peak disappearing 
at Pi + Ap which was present at pi. The resulting momentum distributions are then added 
incoherently, weighed by the height of the momentum distribution of the original BEC at 
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FIG. 3: A sketch of the kick laser setup. A Unearly polarized laser beam is retro-reflected to form 
a standing wave potential. The laser is stabilized to the F = 2 ^ F' = 3 transition in ^^Rb. 

Ap. The final momentum distribution is obtained by the sum of horizontal profiles in Fig. [TT] 
weighed by the initial momentum distribution. 



III. EXPERIMENTAL OBSERVATION OF QUANTUM RESONANCES 
A. Experimental Setup 

We use a double magneto-optical trap (MOT) configuration for the formation of all- 
optical Bose-Einstein Condensate (BEC) of {F = 1) ^''Rb at a temperature of 50 nK. The 
all-optical BEC of ~ 2.10'* atoms is formed in a cross dipole trap created by a pair of 
intersecting focused CO2 laser beams. A detailed description of the experimental setup can 
be found in 
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We realise the AOKR by pulsing a near resonant optical standing wave, derived from a 
780 nm diode laser, onto a BEC of ^''Rb atoms. The kick laser is locked to the Si,F = 2-^ 

2 

Ps, F' = 3 transition in the ^^Rb isotope. Hence, the laser frequency is detuned by 2.45 GHz 
from the relevant F = 1 — t- F = 2 resonance frequency. The laser beam from the diode laser 
passes through an acousto-optic modulator (AOM) for fast switching. After passing through 
the AOM, the beam passes through a single mode optical fibre and is focused onto the BEC. 
The beam diameter at the focus is ~ 100 /im, much larger than the size of the BEC (~10 
/im), thereby yielding a constant interaction strength over the BEC, while increasing the 
intensity. After passing through the center of the trap the linearly polarized beam is then 
retro-reflected to produce a standing wave. The kick laser setup is shown in Fig. 121 

In the AOKR experiment, the BEC is illuminated with short pulses of light in the form 
of a standing wave. The standing wave potential or the interaction potential acts like a 
diffraction grating, which changes the atomic momentum, and the BEC eventually splits into 
a number of momentum components. An alternate picture is that the atoms absorb a photon 
of momentum hk from one beam, and then provide another photon through stimulated 
emission in the other beam. The diffracted components of momentum are therefore separated 
by ±2hkL. In order to reduce mean field effects, the trap that contains the BEC is turned 
off 500 fis before the kick sequence is applied. A shutter is used in the laser beam to ensure 
that it is totally extinguished during the evaporation phase to create the BEC. 

The periodic potential is related directly to the intensity of the laser beam and inversely 
to its detuning from the atomic resonance. High laser power allows the experiment to be 
done at large detunings to suppress spontaneous emission. Therefore, for the interaction 
potential it is desirable to get as much power as possible from the kick laser. The total laser 
power of the experiment can be varied up to 2 mW. The laser pulses are so short that the 
evolution of the wave packet due to its momentum can be ignored during the kicking time. 



35|. 



Atoms satisfying this condition are said to be in the Raman-Nath regime [27|, 

The momentum distribution of the atoms after the kick sequence is measured by absorp- 
tion imaging using the time-of-flight technique, with a typical flight time of 5 to 10 ms. Just 
prior to imaging, the atoms are optically pumped to the F = 2 state by a 100 fis pulse on 
the F = 1 — > F' = 2 repump transition. An absorption image is then obtained using a 
probe laser which is tuned to the 5*1 , F = 2 — Pa , F' = 3 transition. The two dimensional 

2 2 

momentum distributions obtained are summed over the width of the cloud to obtain a one 
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FIG. 4: The momentum distribution of the atoms fitted by a Gaussian distribution to obtain the 
energy. 

dimensional momentum distribution. The energy of the atoms is then obtained by the vari- 
ance of the momentum distribution by fitting a number of Gaussians, one to each diffraction 
order. A typical example of the fitted momentum distributions is shown in Fig. |H 



B. Results of the primary quantum resonances 

The kinetic energy of the system in the case of a standing wave potential (/3 = 0) as a 
function of pulse period is studied for a range of pulses or kicks. This situation can also be 
thought of as equivalent to that of varying the spatial distance between the gratings in the 
optical Talbot effect. To start with, the energy is measured by scanning the pulse period 
(time delay) for two kicks. The period between the kicks is varied from (10 — 70) /is. The 
second kick is used as a tool to study the time evolution of the wave function after the first 
kick. It is well known that in the atomic version of the kicked rotor, at the Talbot time 
Tt = 66.3 fis of free evolution, an identical wave function to that directly after the first kick 
is observed. The introduction of another grating or kick at the Talbot time then simply 
doubles the effect or duration of the first kick. As a result, a quadratic growth in energy is 



observed 
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FIG. 5: Experimental data: (Top) Absorption images for the quantum resonance case (a) and the 
quantum anti-resonance case (b). (Bottom) Energy as a function of the number of kicks (1 to 6) 
for the quantum resonance case (circles) and anti-resonance case (diamonds). The kick strength 
(j)d ^ 2.5 



At a pulse period of ~ 33.1 fis, which is equivalent to ha 
the quantum phase factor between kicks alternates sign 



f of the Talbot time, i.e. Tt/2, 
. As a result, a re- image of the 



grating's transmission function, which is the inverse of the initial one, is obtained. In effect, 
two kicks with half the Talbot time in between them results in the complete negation of the 
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Pulse period (/is) 

FIG. 6: (Color online) The energy as a function of pulse period for two kicks. Blue circles show 
experimental data, whereas the solid green line represents the simulations. The parameters are 
AT = 2, T = (10 - 70) fis, (pd ~ 2.5. 

effect of the first kick. The quantum resonance and anti-resonance behavior is observed in 
the experiment as shown in absorption images in Fig.|5]a and b (Top) with the corresponding 
energy of the atoms representing these resonances (Bottom). 

C. Intermediate secondary resonances 

The energy dependence of the condensate diffracted by two gratings or kicks as a func- 
tion of pulse period between the kicks is shown in Fig. |6l The solid curve represents the 
theoretical simulations, whereas the circles are the experimental observations. The error 
bars are determined by running the experiment a number of times. The well defined peaks 
and valleys that are observed correspond to the integer and fractional Talbot times. The 
maximum in energy at ~ 66 /is is in close resemblance with the simulation, and is a clear 
indication of the Talbot time Tt that corresponds to the quantum resonance. The minimum 
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in energy at ~ 33 /is is also clearly visible, matching the theoretical predictions of half the 
Talbot time Tr/2. 

The energy of the system is also studied as a function of pulse period, by adding another 
pulse (kick) to the first two. The addition of this pulse increases the energy of the system 
by inducing an extra phase shift. The global minimum in this case is then replaced by a 
local maximum as the complete cancelation of the first two kicks is destroyed by the third 
one. This is true for an odd number of kicks. The final energy of the system for three kicks 
as a function of pulse period is shown in Fig. [7] (Top). Most of the structure in Fig. [7] is less 
visible; however, the maximum in energy at about Tt is still obvious as the gratings for this 
amount of free evolution do not induce any additional phase shift. 

The evolution of the system in the presence of large number of kicks shows interesting 
results in the simulations. By measuring the energy, the evolution is examined in the presence 
of four and five kicks. It is well known that with an increase in the number of kicks, the width 
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m. 



of the quantum resonances when plotted against the pulse period becomes narrow 
In Fig. [7] (bottom), the energy as a function of pulse period is plotted for four kicks. In 
this case, at a pulse period of ~ Tt/2, the minimum in energy (close to zero) is observed 
again. As for an even number of kicks, the effect of the previous kick is canceled again for 
the corresponding time delay of Tt/2, and almost no diffraction occurs. The maximum in 
energy at Tt is again visible, since for this time delay the energy of the system increases 
with the addition of further kicks. It should be noted that the resonance peak observed 
at Tt is narrow compared to the peak for a small number of kicks. This is in accordance 
with theoretical simulations, where a narrowing of the quantum resonance peaks for a large 
number of kicks is observed. 

A similar behavior is observed when an extra pulse is added. In Fig. |8l the energy is 
examined by a complete scan of the pulse period up to the second primary resonance for 
five kicks. The minimum in energy at Tt/2 , which was obvious in the case of four gratings, 
however, is replaced by a local maximum as shown. This is because of the extra phase shift 
induced by the fifth grating, due to which the energy increases by some amount for the 
corresponding time delay. It is shown that much narrower resonance peaks are observed in 
this case around the quantum resonances of Tt and 2Tt (~ 132 /is). The maximum in the 
experiment in this case, however, is not in agreement with the simulation. We believe this 
discrepancy is caused by the limited signal to noise ratio. Because of the limited resolution 



14 



40 





20 40 60 

Pulse period (/is) 



80 



FIG. 7: (Color online) The energy as a function of pulse period for three kicks (Top) and four kicks 
(Bottom). The blue circles represents experimental data, whereas the dashed lines are results of 
the simulations. The parameters are T = (10 — 70) /xs, (t)d w 2.5. 
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FIG. 8: (Color online) The energy in recoils as a function of pulse period. Experimental data (blue 
circles), whereas the dashed lines are results of the theoretical simulations. The parameters are 
N = 5, (pd^3 andT = (10 - 140) ^s. 

of the experiment, the observation of the probability of the atoms at high momenta and 
detecting a small amount of diffraction orders at high momenta is difficult, which results in 
the reduction of the energy. Because of finite width of the initial momentum distribution 
the energy of the momentum distribution tends to be extremely sensitive to the population 



at high momenta [36|. Overall, the experimental results are in good agreement with the 
theoretical predictions and simulations. Higher order resonances are expected to occur in 
these systems at fractional multiples of the Talbot time Tt, i.e., at 

m 



T-. 



frac 



-\Tt, 



(18) 



where / and m are integers forming a rational number. The quantum revivals, or resonances, 
at these fractional times are the results of superposition of q copies of the initial wave packet, 
which are separated by 27r/g 4^ . Expecting resonances at these fractional times, the energy 
of the atoms is measured and plotted as a function of pulse period in between the gratings. 
In Fig. [9] the energy is shown for a full scan of the pulse period from (10 — 20) for 5 and 
10 kicks. The maximum in energy at ~ 16.7 /xs is the evidence of a fractional resonance at 
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FIG. 9: (Color online) The energy in recoils as a function of pulse period or time delay between 
the gratings is shown for five (left) and ten (right) kicks with a kick strength of (pd ~ 3. Blue circles 
show experimental data, whereas dashed lines represents the simulations. 

(l/4)Tr. Similarly, the small peak at ~ 13.2 fis shows another resonance which corresponds 
to a fractional Talbot time of {1/5)Tt. 

In Fig. |9] (right), the fractional resonances become sharp as a result of large number 
of kicks, ten kicks in this case. For multiple gratings, the energy of the atoms is more 
sensitive to each induced phase shift and it is extremely difficult to analyze accurately in the 
experiment. The experimental data follow the general trend of those from the simulations. 



IV. INITIAL MOMENTUM DEPENDENCE ON QUANTUM RESONANCES 

In this section, measurements of the effect of the initial momentum on the primary 
quantum resonances in the delta-kicked rotor are discussed. The study of quantum transport 



in such systems has been of great interest in recent years in experiments with cold atoms 41 



421], but its dependence on the initial velocity of the atoms still needs further investigation. 



A. Experimental setup 

In order to take into account the initial momentum dependence on the quantum resonant 
effects, a modified kick laser setup is used, shown in Fig. [TOl The laser beam is split into 
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FIG. 10: A sketch of the laser setup for a moving standing wave: The AFG (Arbitrary Function 
Generator) is used to drive the AOMs, enabhng us to change the frequency of one of the laser 
beam by a small amount. 

two by a beam splitter, and the two beams are then passed through separate AOMs. The 
AOMs used for fast switching are driven in this case by an Arbitrary Function Generator 
(Tektronix AFG-3252), amphfied by home-built RF amplifiers. The AOMs are switched on 
simultaneously, with a tunable frequency difference 5uj. The AFG generates two 80 MHz 
(plus a desired offset frequency), 3.5 Vpp sine waves that shift the frequency of the laser 
beams. The frequency difference gives us an effective initial momentum pi for the atoms, 
given by pi/prec = / {^'^rec)- In the experiment, using the definition of E^ec = h'^k'^/2m, 
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for A = 780 nm one recoil frequency is given by free = 3.77 kHz. From the definition above 
and in [14], 3.77 kHz corresponds to a value of the quasimomentum /3 = .25 recoils, and 
30.16 kHz corresponds to /3 = 2 recoils. An equation editor ArbExpress was used to generate 
a pulse sequence for different pulse periods. The frequency of one beam was set at 80 MHz 
and the second beam was varied by 80 MHz+3.77 kHz to 80 MHz+30.16 kHz in eight steps. 
By introducing the small frequency difference, a moving wave rather than a standing wave 
was obtained. The frequency difference corresponds to quasimomentum values of zero to two 
photon recoils. It should be noted that what is important is the relative phase of the sine 
waves driving the AOMs in the subsequent pulses, caused by the small frequency difference, 
not the actual frequency difference. 

An ensemble of cold atoms with a momentum spread much less than a single photon 
recoil is used. These cold atoms are then kicked at kick periods corresponding to quantum 
resonances. The analytical expression for the amplitudes of the momentum states with 
momentum 2jhk after n kicks has been derived in [3| and [l^, for the kick periods that 
are half integer times the Talbot time Tt, 

= U"^^) Ver^^^^^er--^\ (19) 

where /3 is the quasimomentum in units of 2hkL, and / is an integer, given by / = 2T /Tx- The 
parameter /3 can be used to assign different initial momenta to the atoms in the experiment. 
From equation 16 and 19, the second order moment of the momentum distribution can be 
found as 

oo 

{p') = {2hkLf J2 |c,f(j+/3)'. (20) 

j=-oo 

The values obtained from Eqn. ( fT9l) are found to be in excellent agreement with those 
obtained from the simulation, and with the experimental data. 

B. Results and discussions 

In Fig. [121 absorption images of the atoms are shown for the antiresonance T = Tt/2 case 
after two kicks for different initial momenta. The atoms are allowed to expand for about 
5 ms before they are imaged by the camera, using the time-of-flight method. The different 
initial momentum is imparted to the atoms by setting a small frequency difference between 
the two laser beams. At p = the effect of the first kick is exactly canceled by the second 
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FIG. 11: The probability of finding atoms at the final momenta (horizontal axis) for a range of 
initial momenta (vertical axis) in the simulation. The probability for kick period of T = Tt/2 
(Top) and T = Tt (Bottom) is shown for four kicks, with average initial momentum {pi) — Prec 
and (/irf ~ 1. The original BEC distribution is also shown (as solid curve). 

for the kick period corresponding to half the Talbot time, and we see almost no diffraction, 
as expected from the Talbot effect. At the initial momentum of pi = prec, we see strong 
diffraction, as Tt/2 corresponds to a resonance for these atoms. A similar effect is observed 
for atoms with an initial momentum of pi = Sprec- For atoms with an initial momentum 
of Pi = 2prec and Pi = 'iprec there is again no diffraction. The reason for the absence of 
diffraction in this case is due to the fact that Tx/2 corresponds to an antiresonance for these 
atoms. For the antiresonance case here, the free evolution part of the Floquet operator 



20 





FIG. 12: Absorption images of the atoms after two kicks, with a kick period of T = Tr/2 for 
different initial momenta from to 4 recoils. 



^i(2Ar+i)«5±_Tj,/2 



(21) 



evolve according to 

U free 6 

and depends on the initial momentum. For T7-/2 the free evolution operator depends on 
hk, resulting in resonances and antiresonaces for specific values as described. The laser 
field induces coupling only between momentum eigenstates differing in momenta by integer 
multiples of hk [14!|. 

The dynamics of the system are examined for a range of initial momentum at pulse periods 
of Tr/2 and Tt between the kicks. In Fig. [13] on the left, the momentum distribution of 
the atoms after the kick sequence is examined for the initial momentum pi in 8 steps, from 
to 2 Prec- The resulting ID momentum distributions, averaged over the three repeats 
of the experiment, are shown as a function of the final momentum for two kicks. It can 
be seen that at zero velocity, there is indeed an anti-resonance and very little diffraction 
occurs. For small increments in the initial momentum, small peaks start to appear at higher 
diffraction orders. At Pi = O.bprec, the first and second order diffraction peaks are observed. 
At an initial momentum of Pi = Iprec, the evolution turns into a resonance, with significant 
diffraction. The probability of atoms at higher momenta starts to decrease for pi > Iprec 
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FIG. 13: (Color online) The number of atoms as a function of the final momentum, for a range 
of initial momenta to 2hk. On the left we show the experimental results for two kicks with a 
period of Tt/2. On the right the results are shown for two kicks with a period of Tt- The solid 
red curves represents the experimental data, whereas the simulations are shown as dashed lines, 
for a kick strength of ~ 1- 

and at pi = 2prec the system returns back to an anti-resonance, with very httle diffraction. 
Also shown in the figure are the results of the simulation. When experiments are performed 
with a broad initial momentum distribution, an average over many different Vi yields an 



observed resonance aX T = Tt/I, as has been shown in many publications, see e.g. 



On the right in Fig. [T31 the resulting ID momentum distributions in the experiment and 
the simulations are shown for a period between the kicks of T = T^. The initial momentum 
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FIG. 14: The final energy in recoils as a function of initial momentum from to 2hk for T = Tj'/2 
(left) and T = Tt (right) is shown for the experiment with uncertainties (as markers) and the 
simulation (as solid line) for two kicks, with 0^ w 1. 

Pi is again varied from to 2 p^ec in 8 steps. At an initial momentum of Pi = 0, there is a 
resonance and strong diffraction occurs as expected. The evolution of the system, however, 
turns into an antiresonance at pi = 0.5prec with small diffraction, and back to a resonance 
at Pi = Iprec- The probability of atoms occupying the higher momentum states start to 
decrease again, and at Pi = l.Sp^ec an antiresonance is observed. The evolution finally turns 
into a resonance again at Pi = 2pj.ec, where a strong diffraction is seen. Again, averaging 
over a range of initial velocities will show the same resonance at T = T^- , even though the 
dependence on the initial velocity cycles at twice the rate. At T = 3Tt/2, the cycle of the 
amount of diffraction varying with the initial velocity is at three times the rate we see for 
T = Tt/2, and so on. The simulated curves show good agreement in terms of the relative 
heights of all the diffraction peaks. 

C. Effect of the initial momentum on the kinetic energy 

As discussed earlier, the AOKR dynamics can be examined by the observation of the 
energy or the variance of the momentum distribution in units of the recoil energy Erec = ^^r- 
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FIG. 15: The final energy in recoils as a function of initial momentum from to 2fik for T = 3Tt/2 
is shown for the experiment with uncertainties (as markers) and the simulation (as solid line) for 
two kicks. The kick strength (f)^ ~ 0.9. 

In this section, we study how the energy of the system is influenced by taking into account 
the initial momentum of the atoms. The energy is scanned as a function of the initial 
momentum pi for different pulse periods for a number of kicks. In Fig. [TU the variation of 
the energy on the initial momentum for T = Tt/2 (left) and T = Tt (right) in both the 
experiment and simulation for two kicks is shown. For T = Tt/2 we see a strong peak in 
the energy at about pi = Iprec- Similarly, for the resonant case T = Tt oscillation in the 
energy is observed. These results are consistent with the results presented in Fig. [131 For 
T = 3Tt/2, the energy varies at three times the rate for T = Tt/2, as shown in Fig. [T5l 
For the data points, the averaged profile is fitted to a number of Gaussians, one for each 
diffraction order, to obtain the energies after the kick sequence in each run. The standard 
deviation is then taken to estimate uncertainties in the results. A direct numerical variance 
of the velocity distribution was found to give similar results. Because of the sensitivity to 
small noise peaks at large momenta, these results were less consistent, and were therefore 
not used. 

Next, the energy is examined as a function of the initial velocity for four kicks. As found. 
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FIG. 16: (Color online) The final energy in recoils as a function of initial momentum from to 
2hk for T = Tt/2 (left) and T = Tt (right) is shown for the experiment with uncertainties (as red 
markers) and the simulation (as solid line) for four kicks, with </>^ ~ 1. 

for a large number of kicks, the widths of quantum resonances when plotted as a function 
of kick period get smaller as the number of kicks get larger. Fig. [16] shows the energy 
as a function of the initial velocity for four kicks, for both the anti-resonance T = Tt/2 
(left) and the resonance case T = Tj- (right). On the left, for T = there is an anti- 

resonance at zero momentum, yielding very small energies. There is a small maximum at 
an initial momentum of 1/4 recoils, decreasing again close to zero at 1/2 recoils. At an 
initial momentum of one recoil, there is a stronger maximum in the energy, decaying again 
to small energies at 3/2 recoils before another smaller maximum. At initial momenta of two 
recoils, very small diffraction occurs and the energy returns to a small value close to zero. It 
should be noted that all these features, prominent in calculations, can be reproduced in the 
experiment. For the kick period corresponding to T = Tj-, a larger number of oscillations 
in energy in the simulation is observed. The maximum in these curves is observed to be 
much narrower than those for a small number of kicks. The maxima in the experiment is 
in agreement with that in the simulation, but the small-period oscillations in the simulation 
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FIG. 17: (Color online) The momentum distributions of the atoms after four kicks. The momentum 
distributions for the kick period of T = Tx/'i (Top) and T = Tt (bottom) are shown in the 
experiment (solid blue lines) and the simulation (red dashed lines) for T = Tt/2. The initial 
momentum pi = prec in both cases and (j)d ^ 1- 

were not resolved in the experiment due to limited resolution for these settings. 
D. Details of the momentum distribution 

In a similar experiment, the momentum distributions obtained as a result of four kicks 
are shown in Fig. [17] for the antiresonance case (top) and the resonance case (bottom). A 
change in widths of the momentum distributions for different diffraction orders is observed 
in both cases. The observed difference in widths can easily be related to the separation of 
one velocity peak into multiple peaks in the simulation. For instance, the velocity peak at 
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a momentum p = 1 recoil is split into three, with varying heights in the two situations. 
This would be caused by the peak disappearing for Ap = 0.2 recoils, and subsequently re- 
appearing for Ap = 0.4. In the final momentum distribution, the latter peaks appear even 
though they are not strongly weighted by the initial velocity distribution. 

Similarly, the peaks at p = 3 and p = —1 recoils are split into two. This would be caused 
by the peaks appearing at a small Ap. In addition, the peaks at p = —(5, 7) recoils are found 
to be considerably narrower than the initial momentum distribution, which is caused by the 
fact that this diffraction peak only appears for very small Ap. These observations can be 
verified in Fig. [TT]by drawing a horizontal line at fixed Ap and finding the diffraction orders. 
In the simulation, it has been found that the width of resonances as a function of the initial 
momentum gets smaller in the presence of a large number of kicks. 

V. CONCLUSIONS 

In this paper we present an experimental study of quantum resonances in the 5— kicked 
rotor system in detail. Quantum resonances are intrinsic to quantum mechanics and are 
therefore useful in studies related to the phenomenon of quantum chaos. We believe that 
the results provided in this paper contribute to a deeper understanding of quantum dy- 
namics. Using the kinetic energy measurements, quantum resonances and antiresonances 
have been observed for a range of kicks in an AOKR experiment. In particular, we have 
shown how these resonances behave as we increase the number of kicks. For large number 
of kicks, the energy of the atom is very sensitive to each induced phase shift and is diffi- 
cult to analyze accurately because of limited experimental resolution. Within the limits of 
experimental resolution we have found that these resonances get sharper as the number of 
kicks is increased, in good agreement with simulations. Moreover, higher order resonances 
that arise from fractional quantum revivals have been observed for pulse periods of {1/4)Tt 
and (l/5)Tr for five and ten kicks. 

The study of quantum resonant effects in this work has been further extended to analyze 
the effect of the initial momentum of the atoms on quantum resonances in the AOKR. The 
first three quantum resonances have been examined for two and four kicks. A sinusoidal 
dependence of the energy on the initial momentum for two kicks has been observed. By 
increasing the number of kicks, the system becomes extremely sensitive to the initial mo- 
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mentum of the atoms. A more complex structure of the energy is observed for four kicks. 
For four kicks, the maxima in the experiment is in agreement with that in the simulation. 
The small-period oscillations in the simulations, however, were not resolved in the exper- 
iment due to limited resolution for these settings. The maximum in the energy curves is 
observed to be much narrower than those for small number of kicks, which is again in good 
agreement with simulations. With these experiments, it has been found that by applying a 
small frequency difference between the beams constituting the standing wave, we can dial 
any initial velocity we choose for the atoms with respect to the standing wave. In general, 
we found that the energy is periodic with initial momentum. 

In the future it would be an interesting study to examine the dependence of the initial 
velocity of the atoms on the fractional resonant effects in the AOKR. In the "ballistic" 
regime, i.e. for the kick period close to one of the "quantum resonances" , a recent proposal 



17( 1 discusses the re-creation of the wave function by a strong, final pulse. The aim would 
be to confirm their predictions and investigate the influence of decoherence. In particular, it 
would be interesting to investigate the influence of the non-linearity caused by having a dense 
atom cloud. By manipulating the relative phase of the delta-kicks, a primitive "quantum 
computer" was recently demonstrated in [43]. The aim would be to take this idea further, 
making use of the excellent control we have of the kick laser field. As these systems are 
very sensitive to initial conditions, the aim would be to investigate and measure these initial 
conditions, in particular to set an initial velocity and measure the recoil frequency, which 
is proportional to where m is the mass of the particle. This was outlined in a recent 



study 144 1 . 
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